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In this paper, we will discuss the prospects of deducing the proton-air cross section from fluorescence telescope 
measurements of extensive air showers. As it is not possible to observe the point of first interaction X\ directly, 
other observables closely linked to X 1 must be inferred from the longitudinal profiles. This introduces a dependence 
on the models used to describe the shower development. The most straightforward candidate for a good correlation 
to Xi is the depth of shower maximum X max - We will discuss the sensitivity of an X max -based analysis on (J p _ a i r 
and quantify the systematic uncertainties arising from the model dependence, parameters of the reconstruction 
method itself and a possible non-proton contamination of the selected shower sample. 



1. Introduction 

One of the most important parameter needed 
for interpretation of air shower measurements is 
the proton-air cross section cr p _ a i r , which is not 
well known at cosmic ray energies. Collider mea- 
surements of the proton- (anti)proton cross sec- 
tion have to be extrapolated over a wide energy 
range to be used in models of nuclear effects to 
predict the proton-air cross section at ultra high 
energies. On the other hand, it might be possi- 
ble to measure er p _ a ir by analyzing extensive air 
showers (EAS). 

The point of first interaction in slant depth X\ 
follows the exponential distribution 

dP ^) 1 e -*/W (i) 

dXl Ap_ a i r 

The cross section cr p _ a i r is related to the absorp- 
tion length A p _ air as A p _ air = M mcan /cr p _ air , with 
M mca ,n being the mean mass of a target (air) par- 
ticle. 

The depth and characteristics of the first inter- 
action of the primary cosmic ray particle with 
a molecule of the atmosphere determines to a 
large extent the resulting shower shape. This im- 
plies strong correlations of the first interaction 
with experimentally observable parameters of the 
shower. 

Several air shower experiments have provided 
measured values of <7 p _ a i r by exploiting correla- 



tions of X\ with observables found by using air 
shower simulations }l|2|3|4|5j . Of course, these 
results depend strongly on the used simulations 
and especially on the chosen high energy hadronic 
interaction model [617] . 

The HiRes collaboration proposed a new tech- 
nique which appeared to be almost independent 
of the hadronic interaction model [5] . The idea is 
to decompose the measured distribution of shower 
maxima <iP(X max )/<iX max into the two indepen- 
dent distributions of X\ and AX with 

X max = X 1 + AX. (2) 

Using the introduced AX one can express the 
X max distribution as folding integral 

poo 

P(X max ) = / dXi Ppfi) • P(AX\X!) (3) 
Jo 

of the Xi and AX distributions, whereas in the 
HiRes [8, approach P(AX\Xi) = P(AX) is used. 

Measuring er p - a ir means to deconvolve the 
Xmax distribution using a known AX distribu- 
tion. 

The distribution of AX needs to be inferred 
from simulations. 

2. Simulation study 

To investigate the deconvolution technique, 
we simulated sets of proton-induced showers of 
10 EeV using CONEX [5] with the high energy 
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Figure 1. Left panel: Dependence of the AX-distribution on Xi. Right panel: Interaction model 
dependence of AJT-distribution. 



interaction models QGSJET01 [ID], SIBYLL2.1 
[II], NEXUS3 [12] , and QGSJETII.3 [13]. In 
addition helium, CNO and gamma-ray primaries 
were simulated with QGSJET01. Each set of sim- 
ulations contains 200.000 events. 



3. Universality of AX distribution 

One of the important conditions for the appli- 
cability of the method proposed in [5] is the X\- 
indcpendence of the AX distribution. For a given 
model, this is indeed the case as shown in Fig. [I] 
(left). The universality in X\ has been verified 
for proton primaries using NEXUS3, SIBYLL2.1, 
QGSJETII.3, and QGSJET01 at 1 EeV and 
10 EeV. 

Therefore one can use 



P{AX\X X ) = P(AX). 



(4) 



However, the resulting AX distribution de- 
pends strongly on the hadronic interaction model 
used for shower simulation (Fig. [I] right). The 
maximum of the AX distribution is shifted by up 
to 60 g/cm 2 and the slope beyond the maximum 
changes by up to a factor of 1.36. This introduces 
a principal model dependence to the deconvolu- 
tion technique. It can be tried to reduce this de- 
pendence by adding an additional degree of free- 
dom, JTshift, which allows the AX distribution to 



shift along the x axis 
P(AX) ^ P(AX + X shm ). 



(5) 



It is much more difficult to reduce the model- 
dependence of the shape of the AX distribution. 
This will not be discussed in this paper but we 
shall show the impact of this model-dependence 
on the resulting £J p _ a i r reconstruction. 

4. Study of the deconvolution method 

For each model, a Monte Carlo data set of 
50.000 events was generated and analyzed with 
the AX distributions found previously. No de- 
tector effects are considered. Two free parameters 
(Ap-air and X^ift) are fitted using x 2 minimiza- 
tion and 



dP(X max ) 



dXi 



1 



■Xa.min ^p— air 

f dP(AX + X shih ) 
\ dAX 



■ (6) 

model 

The fit results are shown in Tab. Q] for all com- 
binations of hadronic interaction models and AX 
distributions. The diagonal elements of the ta- 
ble represent matching models for data and AX. 
In these cases the true cross section of the input 
model can be retrieved very well by the convo- 
lution method. It is interesting that for mod- 
els with a similar shape of the AX distribution 
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Table 1 



Analysis results of cr p _ a i r in mb for 10 EeV proton showers. The results for Y sn ift are given in g/cm 2 . 
The expected cross sections are 532mb (QGSJETOl), 612mb (SIBYLL 2.1), 542mb (NEXUS), and 
562 mb (QGSJETII). 



data 


A X model used for reconstruction 


data model 


QGSJET 

free ^shift 


SIBYLL 

Croc ^shift 


NEXUS 

free ^shift 


QGSJETII 

free ^shift 


QGSJETOl 
SIBYLL 
NEXUS 
QGSJETII 


537.1±5.5 1.7±0.8 
802.8±12 36.1±0.7 
749.9±10 55.1±1.2 

697.1±8.2 27.0±0.8 


465.3±4.1 -29.8±0.8 
609.2±5.7 0.7±0.5 
569.1±3.9 19.5±0.9 
553.2±4.3 -7.5±0.6 


447.9±3.2 -46.9±0.8 
572.8±4.5 -17.4±0.5 
543.1±4.0 1.5±0.5 
521.6±3.7 -26.1±0.3 


467.6±3.7 -19.6±0.7 
613.3±5.7 10.6±0.5 
576.5±4.6 29.5±0.8 
562.3±4.4 2.1±0.4 



(QGSJETII and SIBYLL) the reconstruction still 
works rather well, whereas for combinations of 
models which have a very different shape of AY 
(QGSJET and NEXUS) the results differ up to 
200 mb. 

4.1. Importance of fit range 

The reconstructed value of <7 p - a i r is very sensi- 
tive to the chosen starting point of the fit range, 
as shown in Fig. [2] (left). The reason for this is 
the introduced Y sn ift parameter which needs the 
peak of the Y max distribution to be restricted. 
To fit only the slope after the peak of the Y max 
distribution is not sufficient if Y sn ift is treated as 
free parameter. 

Restricting the end point of the fitting range 
mainly increases the statistical uncertainty of the 
fit result. The true value of er p ™| lr is still within 
the statistical errors of the reconstructed values 
(see Fig. [2j right). 

4.2. Impact of other elements 

If the selected showers do contain a fraction of 
light elements other than protons, the shape of 
the measured Y max distribution will be altered. 
A possible contribution of helium will mostly dis- 
tort the region of the peak of the Y max distribu- 
tion. Since the peak of the Y max distribution is 
important for the fit itself, the reconstruction of 
(Tp-air is affected (Fig. O top) 

Contamination with CNO showers mainly in- 
fluences the rising edge of the Y max distribution. 
The impact on the fit result can be limited by 
shifting the starting point of the fit range more 
towards the peak of the distribution. Using a 
short fit range, the effect on the reconstruction 
of cTp_ a i r can be neglected for small fractions of 
CNO (Fig. El middle). 



4.3. Impact of gamma-ray primaries 

Gamma-ray primaries produce an Y max distri- 
bution which is shifted to higher Y max by about 
200 g/cm 2 with respect to protons. Therefore the 
influence of gamma-ray showers can be reduced 
by setting the end of the fitting range appropri- 
ately (Fig. [3l bottom). Nevertheless already a 
small contribution of gamma-rays does affect the 
Y max distribution significantly because of the low 
number of proton events at high Y max . It should 
be noted that even if cosmic ray sources would ac- 
celerate only protons, a small fraction of gamma- 
rays is expected due to the GZK energy loss pro- 
cess [14] . 

5. Summary 

The folding integral method would be very 
powerful given the true AY distribution is 
known. Input cross sections can be reconstructed 
within statistical uncertainties of a few mb. Un- 
fortunately, the AY distribution does depend on 
the high energy hadronic interaction model used 
during air shower simulation. The AY model de- 
pendence can be reduced by adding a new degree 
of freedom Y sn ift to the fit but it still amounts 
up to ~ 100 mb. On the other hand, treating 
Y s hift as free parameter increases the sensitiv- 
ity of C7p- a i r to a possible helium contamination 
of the Y max distribution. Medium mass cosmic 
ray primaries have only limited influence on the 
reconstructed £r p _ a ir. Even a small fraction of 
gamma-ray primaries has a significant impact on 
the reconstruction results. 

To use the folding integral method for a reli- 
able cross section measurement, the model depen- 
dence of the AY distribution has to be reduced 
significantly. 
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Figure 2. Impact of the begin (left panel) and the end (right panel) of the fit range of the X max distribution 
(pure proton composition at 10 EeV). On the x axis the distance to the peak of the A max distribution is 
shown. The horizontal lines denote the true value of Op™°j r modcl for each interaction model. 
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Figure 3. Impact of non-proton contributions in 
the analyzed shower sets on the reconstruction of 
(7p_ a i,.. The horizontal lines denote the true value 

r true, model r , i i 1 • i ■ j_- 

ot er p _ a j r ot the used mgn energy interaction 
model. 



